Elliptic flow of thermal photons in Au+Au collisions at ^/sNN = 200 GeV 
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The transverse momentum (pt) dependence, the centrahty dependence and the rapidity depen- 
dence of the elliptic flow of thermal photons in Au+Au collisions at ^/snn = 200 GeV are predicted, 
based on a three-dimensional ideal hydrodynamic description of the hot and dense matter. The el- 
liptic flow parameter V2, i.e. the second Fourier coefficient of azimuthal distribution, of thermal 
photons, first increases with pt and then decreases for pt > 2 GeV/c, due to the weak transverse 
fiow at the early stage. The pt-integrated V2 first increases with centrality, reaches a maximum at 
about 50% centrality, and decreases. The rapidity dependence of the elliptic fiow V2{y) of direct 
photons (mainly thermal photons) is very sensitive to the initial energy density distribution along 
longitudinal direction, which provides a useful tool to extract the realistic initial condition from 
measurements. 

PACS numbers: 



I. INTRODUCTION 

The deconfined and novel nuclear matter, the quark 
gluon plasma (QGP), has been expected to appear in 
relativistic heavy ion coUisions. Various signatures have 
been proposed to verify its existence Many experi- 
ments have been done so far to explore properties of the 
QGP at the Relativistic Heavy Ion Collider (RHIC) and 
will be done at the Large Hadron CoUider (LHC). Ex- 
perimental results indicate that the QGP is created at 
RHIC 13. 

Collective flow, in particular eUiptic flow, is an effec- 
tive probe to investigate bulk properties of the QGP in 
nucleus-nucleus collisions at relativistic energies. In non- 
central collisions the overlapping reaction zone of two col- 
liding nuclei has an anisotropic shape, like an almond, in 
the transverse plane. This leads to a preferred flow direc- 
tion parallel to the impact parameter, and therefore the 
initial spatial anisotropy is carried over to momentum- 
space anisotropy The QGP created at the early stage 
expands, cools down, and flnally forms a hadronic gas 
(HG) phase. Then hadrons still have strong rescatterings 
until freeze-out. So those hadrons do carry information 
on the asymmetric flow, but only concerning the freeze 
out surface. 

Fortunately, in high energy heavy ion coUisions, ther- 
mal photons can be produced during the whole history 
of the evolution of the hot and dense matter. Moreover 
the mean free path of photons is much larger than the 
transverse size of the bulk matter. So thermal photons 
produced in the interior of the plasma pass through the 
surrounding matter without any interaction. As a re- 
sult, thermal photons provide information on flow asym- 
metries even from inside the bulk volume, not just its 
surface. 

In this paper, we investigate the relation between the 
observables of thermal photons and the space-time evo- 
lution of the hot and dense matter for exploring QGP 



properties, based on a fully three dimensional (3D) ideal 
hydrodynamic description of expansion of the matter. 
Toward establishment of the relation, we calculate the 
second Fourier coeflicient of azimuthal momentum distri- 
bution, the so-called eUiptic flow coefficient D2, for ther- 
mal photons in Au+Au collisions at ^/snn — 200 GeV 
at various centralities (0-70%). 

In Sec. 2 we will briefly review the space-time evolu- 
tion of the hot and dense matter using the 3D ideal hy- 
drodynamics and the basic formula for the production of 
thermal photons. In Sec. 3, we will show our results on V2 
of thermal photons in Au+Au collisions at ^/sNN = 200 
GeV, its transverse momentum, centrality dependences 
and rapidity dependence. Section 4 is devoted to discus- 
sion and summary of our results. 



II. BULK EVOLUTION, THERMAL PHOTONS, 
AND ELLIPTIC FLOW 



In our calculation, a full 3D ideal hydrodynamic calcu- 
lation [1, Ol is employed to describe the space-time evolu- 
tion of the hot and dense matter created in Au+Au colli- 
sions at RHIC. The local thermal equilibrium is assumed 
to be reached at the initial time tq = 0.6 fm/c. The 
transverse flow is assumed to be zero at tq . We consider 
two scenarios to obtain entropy density (or energy den- 
sity) at TTi: one is a parameterization based on Glauber 
modeljl, @|, and the other is based on flux tubes (string 
picture) reahzed in the EPOS modelf^. For these two 
initial conditions, we take the same hydrodynamic equa- 
tions and equations of state. The initial energy density 
from the two scenarios has been plotted as a function 
space-time rapidity rjs in FigUl where the dashed lines 
refer to the parameterized initial conditions, and dotted 
lines to the EPOS case. One can see that the energy 
density decreases very rapidly with r]s in EPOS, while 
it has a plateau for \r]s\ < 2.5 using parameterized ini- 
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Figure 1: (Color online) energy density as a function of space- 
time 77s ;Dashed lines from parameterized initial condition and 
dotted lines from EPOS initial condition. 



tial conditions. In two dimensional hydrodynamics this 
plateau is assumed for the whole rjg region. The EPOS 
initial condition gives many similar results compared to 
the parameterized ones, i.e., the measured pt-spectra of 
direct photons can also be well reproduced, and the same 
Pt-dependence and centrality dependence of elliptic flow 
of thermal photons has been observed. But there are also 
striking differences, to be discussed later. 

For both scenarios, the impact parameters correspond- 
ing to different centralities in Au+Au collisions at RHIC 
are 3.2, 5.5, 7.2, 8.5, 9.7, 10.8, and 11.7 fm for 0-10%, 
10-20%, • • • and 60-70% centrality, respectively. The 
evolution of the fluid is governed by conservation laws of 
energy and momentum 



^^T^'•' = 



with the energy-momentum tensor 

T"" = (e + P)m^u'' - Pg' 



(1) 



(2) 



for a perfect fluid. Here, e, P, and m'' are energy den- 
sity, pressure, and local four fluid velocity, respectively, 
which are related through the equation of state (EOS). 
The above equations are solved in full 3D space (r, x, y, 
rjs) where r, 77s, x, and y are the proper time, space-time 
rapidity, the two transverse coordinates along the im- 
pact parameter and perpendicular to the reaction plane, 
respectively. Since the net-baryon density is small near 
the mid-rapidity at RHIC, it is neglected. We also ne- 
glect the dissipative effects in the space-time evolution. 
The critical temperature of a flrst order phase transition 
between the QGP phase and the hadron phase is flxed at 
Tc = 170 MeV. The EOS for the QGP phase is 



P = -(e-4B), 



(3) 



where the bag constant is tuned to ~ 247.19 MeV 
so that pressure of the QGP phase with three flavors is 
matched to that of HG phase (up to mass of A(1232)) 
at Tc. In the hadronic phase, we employ two resonance 



gas models: partial chemical equilibrium (PCE) and full 
chemical equilibrium (FCE) In the PCE model, the 
particle ratios including contribution from resonance de- 
cays are assumed to keep during the hadronic expan- 
sion, which is required from the data of particle ratios. 
The PCE model is our standard choice to calculate the 
photon spectra in the following. For the purpose of 
comparison, we also employ the conventional resonance 
gas model in which chemical equilibrium is maintained 
during the hadronic expansion. In both models, the 
hadronic matter is kinetically and thermally frozen at 
e'^ = 0.08 GeV/fm'^. Corresponding freezeout tempera- 
ture is T'h - 100 (130) MeV for the PCE (FCE) model, 
which is consistent with the pt slope of protons. 

Transverse momentum spectra of thermal photons can 
be written as 



dN 
dy d^pt 



= j d'^xT{E\T) 



(4) 



with V{E* ,T) being the Lorentz invariant thermal 
photons emission rate which covers the contributions 
from the QGP phase Q and HG phase d^x = 
T dr dx dy drjs being the volume-element, and E* = p^u^ 
the photon energy in the local rest frame. Here p^ is 
the photon's four momentum in the laboratory frame, T 
and are the temperature and the local fluid velocity, 
respectively, to be taken from the hydrodynamic calcula- 
tions mentioned above. We only consider thermal photon 
radiation from the region with energy density above e'^. 
For more details one can check where the measured 
Pt spectra of direct photons at different centralities in 
Au+Au collisions are nicely reproduced. 

The triple differential spectrum can be written as a 
Fourier series. 



d^N 

dy d^pt 2npt dpt dy 



n=l 



2vn cos{n(j)) 



(5) 



where (j) is the azimuthal angle of photon's momentum 
with respect to the reaction plane, which is deflned to 
be the plane containing the impact parameter and beam 
axis. The elliptic flow is quantifled by the second har- 
monic coefficient 1)2 



V2{pt,y) 



J d(t)Cos{2(j))d^N/dyd'^pt 
Jd(j)d^N/dyd^Pt ' 



(6) 



The Pt dependence of the triple differential spectra is 
strongly affected by the flow u through the argument 
E* ~ p'^u^ in the photon emission rate. In the local 
rest frame thermal photons of any given energy are emit- 
ted isotropically. The isotropic distribution is Lorentz- 
boosted by flow velocity zv. The azimuthal asymmetry of 
the transverse components of the flow obviously ends up 
with an anisotropic momentum distribution, which gives 
a finite V2. Therefore both strength and the anisotropy 
of transverse flow velocity are important to generate the 



3 



elliptic flow of thermal photons. To understand the ther- 
mal photon V2 quantitatively, we define the mean radial 

flow {Vr} and the mean anisotropy of fiow (^v^^'^^"'^ as 

M = (yfif^y) (7) 
= = (8) 

where (• • • ) stands for energy-density- weighted space- 
time average, and Vy are the flow velocity components 
along X-axis and y-axis, respectively. We only consider 
the region above e — e'^ = 0.08 GeV/fm'^, where thermal 
photons are produced in this model. 

Table I: The strength and the anisotropy of transverse flow 
at each centrality. 

Centrality(%) 0-10 10-20 20-30 30-40 40-50 50-60 60-70 
(vr) 0.114 0.122 0.123 0.117 0.109 0.0959 0.0804 

^hydroX 0.0417 0.103 0.154 0.188 0.212 0.222 0.240 



From Table HI we see that the mean radial flow in- 
creases with centrality from 0-10% {b = 3.2 fm) to 20- 
30% {b = 7.2 fm), then decreases from 20-30% (5 = 7.2 
fm) to 60-70% (6 = 11.7 fm). On the other hand, the 
mean flow anisotropy increases with centrality from 0- 
10% to 60-70% monotonically. One can easily understand 
the tendency when we investigate the time evolution of 
space-averaged results in Figs. [2] (a) and (b). One should 
notice that, for all centralities, the average radial flow 
is quite small before t = 1 fm/c during which most of 
thermal photons in 4 < < 6 GeV/c are produced. 

III. RESULTS 

In Fig. O the V2 of thermal photons in < pt < 6 
GeV/c is shown for various centralities from to 70% in 
Au+Au collisions at ^/snn = 200 GeV. The solid lines 
from bottom to top refer to centralities 0-10%, 10-20%, 
20-30%, 30-40%, and 40-50%, respectively. The dashed 
lines from top to bottom at pt = 2 GeV/c refer to the 
centralities 50-60% and 60-70%. 

For each centrality, the thermal V2 increases then de- 
creases with increasing pt and a peak appears at pt ~ 
2 GeV/c. The same pt dependence has been predicted 
in the study with 2+ ID ideal hydrodynamics [ll|. The 
decrease of thermal photons' V2 at high pt can be ex- 
plained as the weak transverse flow at the early stage. 
Because at the higher pt, the more fraction of thermal 
photons are emitted from the hot matter at an early stage 
for all centralities. i.e., at pt = 3 GeV/c, about 50% 
thermal photons are produced during r £ (0.6, 0.9)fm/c 
and at Pt = 4 GeV/c, the fraction is about 70%. In 
our 3+ ID hydrodynamics, the evolution of the radial 
flow is plotted in Fig. [2ja). At an early time near tq. 



i.e.,T S (0.6, 0.9)fm/c, the transverse flow vanishes, and 
so does the elliptic flow of thermal photons. Therefore, 
the decrease of thermal photons' V2{pt) at high pt reflects 
the weak transverse radial flow in the early stage. Note 
that a similar non-trivial behavior is also seen in pt slope 
parameters of dimuon spectra as a function of invariant 
mass [lo| . 

Similar to the 2+ ID hydrodynamics(Tl| . small bumps 
appear in the elliptic flow of thermal photons at pt close 
to zero. This is because of the thermal photons are emit- 
ted from two phases, QGP phase and hadronic phase at 
different time. In fact, thermal photon emission rate 
is not reliable at pt — > 0. 

One can also see an EOS dependence of thermal pho- 
ton elliptic flow. As an example, we change from the 
PCE model to the FCE model, which has been employed 
in the conventional hydrodynamic calculations including 
the thermal photon calculations in the earlier study [ll|. 
The change of the EOS is almost invisible in the pt spec- 
tra of thermal photons due to the dominant contribution 
from the QGP phase. However, V2{pt) from the hadronic 
phase is about ten times bigger than the one from the 
QGP phase. So a slight change in the hadronic phase 
can be magnified and become quite visible in the elliptic 
flow of thermal photons, as shown in Fig. |4l One can 
also conclude that the elliptic flow of thermal photons is 
sensitive to the EOS, particularly in the hadron phase. 
However, the effect of different EOS we obtained here is 
quite similar to the effect from the different formation 
time of quark gluon plasma • 

The centrality dependence has been investigated in 
FiglU where (a) represents the yield and (b) the V2 
of thermal photons at midrapidity, and where different 
types of curves present different pt range. It is clear that 
the yield of thermal photons decreases with centrality 
monotonously. However, the elliptic flow of thermal pho- 
tons does not change with centrality monotonously. It 
reaches a maximum at impact parameter b=9.7fm or 40- 
50% centrality, then decreases for more central or more 
peripheral collisions, due to the interplay between the 
strength and the anisotropy of transverse flow. As we 
discussed in Sec. 2, both the strength and the anisotropy 
of transverse flow are important factors to generate the 
elliptic flow of thermal photons. The anisotropy increases 
with decreasing centrality, but the strength of transverse 
flow decreases, as shown in Table [B 

In FigH] the pt ranges used are (0,3)GeV/c: dotted 
lines; (0,6): solid lines; (0.5, 6): dashed lines; (1,6): 
dashed dotted lines. One can see the yield and the pt- 
integrate f 2 of thermal photons are not sensitive to the 
upper limit, c.f. the integral range (0,3) and (0,6) GeV. 
This is clear because the pt-spectrum of thermal photons 
decrease almost exponentially and high pt region con- 
tributes very little to the total yield. The pt-spectrum 
of direct photons is dominant by the thermal photons 
at low Pt region and also decreases rapidly with pt. So 
we can expect that the pt-integrated quantities of direct 
photons behaviors similarly to those of thermal photons. 
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Figure 2: (Color online) Average radial flow velocity (a) and average flow anisotropy (b) as a function of proper time. Different 
curves correspond to 0-10%, 10-20%, ... 60-70%, respectively. A line with longer life time represents the corresponding result 
for a more central collision. 
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Figure 3: (Color online) V2 for thermal photons in Au+Au 
collisions at y^SNN = 200GeV is shown for various central- 
ities from to 70% in < pt < 6 GeV/c. The solid lines 
from bottom to top correspond to centralities 0-10%, 10-20%, 
20-30%, 30-40%, and 40-50%, respectively. The dashed lines 
from top to bottom at pt = 2 GeV/c correspond to the cen- 
tralities 50-60% and 60-70%. 
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Figure 4: (Color online) The elliptic flow of thermal photons 
at & = 7.2 fm. Solid (Dashed) line is the result from a par- 
tial chemical equilibrium model (a full chemical equilibrium 
model). 



On the other side, one can see that the yield and the pt- 
integrate V2 of thermal photons are very sensitive to the 
lower limit, due to the rapidly decreasing pt-spectrum. 
One should be careful when a comparison between the 
prediction and the measurement is performed since the 
lower limit is experimentally determined by the detectors. 

The above results are all at midrapidity. In the fol- 
lowing we discuss the rapidity dependence, for the two 
scenarios: initial conditions based on EPOS flux tubes 
and the parameteraized one. As already said, both sce- 
narios give many similar results. However, the rapidity 
dependence of elliptic flow is very different: in FiglHl the 
rapidity distributions oidn/dy and the V2 of thermal pho- 
tons produced from the two kinds of initial conditions are 
shown. Dashed lines for parameterized one and dotted 
lines for EPOS. One can see thermal photon yield dn/dy 
has a similar dependence on rapidity y, but two very dif- 
ferent dependences of the elliptic flow V2{y) have been 
obtained. One can see that V2{y) shows almost the same 
shape as the initial e(??s): A rapid decrease of eUiptic flow 
along longitudinal direction is obtained from EPOS ini- 
tial conditions, where the energy density decreases very 
rapidly with rjs, while a plateau of V2{y) at midrapid- 
ity region is obtained from the parameterized initial con- 
ditions. A similar observation has already been made 
concerning the rapidity dependence of elliptical flow of 
hadronsQ. 



IV. DISCUSSION AND CONCLUSION 

Based on a fully three-dimensional ideal hydrodynamic 
description of the evolution of hot and dense matter cre- 
ated in Au+Au coUisions at ^Jsnn = 200 GeV, we found 
that the eUiptic flow of thermal photons depends non- 
trivially on not only the anisotropy but also the strength 
of transverse flow. 

At midrapidity, the V2 of thermal photons increases 
with pt up to pt ~ 2 GeV/c, then decreases at higher 
Pt, similar to the prediction from two-dimensional ideal 
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Figure 5: (Color online) (a) the yield and (b) the V2 of thermal photons at midrapidity. Different curve types represent different 
Pt-integral ranges, i.e. (0,3)GeV/c: dotted lines; (0,6): solid lines; (0.5, 6): dashed lines; (1,6): dashed dotted lines. 
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Figure 6: (Color online) (a): rapidity dependence of thermal photon production, (b): rapidity dependence of the pt-integrate 
V2 of thermal photons. Dashed lines from parameterized initial condition and dotted lines from EPOS. 



hydrodynamics . The decrease of thermal photon V2 
at high pt reflects the fact that at very early times, the 
transverse radial flow is weak. So thermal photons can 
serve as a penetrating and unique probe for the whole 
history of the evolution. 

The pt-integrated V2 of thermal photons increases with 
increasing impact parameter 6 up to ~ 10 fm (40-50% 
centrality) and then decreases above. Thermal photons 
from two different phases have the same centrality depen- 
dence. The decrease of thermal photon V2 at peripheral 
collisions is due to the interplay between the strength and 
anisotropy of transverse radial flow of the hot matter: As 
increasing impact parameter b, the anisotropy increases, 
but the strength of transverse radial flow decreases. 

The rapidity dependence of elliptic flow of thermal 



photons V2{y) can "remember" the initial energy density 
along the longitudinal direction eo{r/s), i.e., a plateau 
at midrapidity can be obtained from the parameterized 
initial condition based on Glauber model, while a rapid 
decrease of elliptic flow along rapidity is obtained from 
EPOS initial conditions. 
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